JOURNAL ON SELECTED AREAS IN COMMUNICATIONS 



1 



Diversity Limits of Compact Broadband 
Multi-Antenna Systems 

Pawandeep S. Taluja, Member, IEEE, and Brian L. Hughes, Member, IEEE 



O 

(N 

> 

O 

O: 

m 
(N 



C/3 



> 

On- 

^ " 



X: 



Abstract — In order to support multiple antennas on compact 
wireless devices, transceivers are often designed with matching 
networks that compensate for mutual coupling. Some works have 
suggested that when optimal matching is applied to such a system, 
performance at the center frequency can be improved at the 
expense of an apparent reduction in the system bandwidth. This 
paper addresses the question of how coupling impacts bandwidth 
in the context of circular arrays. It will be shown that mutual 
coupling creates eigen-modes (virtual antennas) with diverse 
frequency responses, using the standard matching techniques. We 
shall also demonstrate how common communications techniques 
such as Diversity-OFDM would need to be optimized in order to 
compensate for these effects. 

Index Terms — multiple antennas, mutual coupling, broadband 
matching, MIMO, OFDM. 



I. Introduction 

MULTIPLE-antenna systems have been shown to alle- 
viate the problem of signal fading, and promise high 
spectral efficiencies in wireless propagation environments rich 
in multipath 0], Q. However, most of these advantages from 
multiple-input multiple-output (MIMO) systems can only be 
realized for large antenna spacings. As the form-factor of 
current handheld and portable devices continues to shrink, and 
the demand for high data rate systems continues to grow, it is 
becoming increasingly necessary to deploy multiple antennas 
in a small space. Several wireless standards, including 4G LTE, 
IEEE 802. 1 In, provide support for multi-antenna devices such 
as mobile phones, laptops, tablets or access points - some of 
which offer a great fit for two dimensional arrays, popularly 
circular. 

With closely spaced antennas, impairments such as fading 
correlation and mutual coupling become increasingly domi- 
nant. Transceivers designed for compact antenna arrays often 
employ special radio-frequency (RF) networks called matching 
networks that are embedded between the antenna array and 
the rest of the RF chain designed for optimized performance. 
Several studies have proposed optimal transceiver design for 
MIMO systems in the presence of mutual coupling by the 
use of multiport matching networks I!]-!!]. However, the 
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focus of these studies has largely been narrowband systems; 
the effects of mutual coupling on the bandwidth of MIMO 
systems has received little attention. These results suggest that 
when optimal matching is applied to a system with strong 
mutual coupling, performance at the center frequency can be 
improved at the expense of an apparent reduction in the system 
bandwidth [9|, ifTOl . This raises fundamental questions about 
the physical realizability of these networks and the bandwidth 
assumptions. 

The study of matching networks for wideband systems 
falls under the purview of broadband matching theory - a 
field rich in multiport broadband matching network design, 
albeit, limited to systems not involving coupled sources or 
loads. However, there exist techniques in antenna/microwave 
circuit design that deal with coupled systems by use of 
orthogonal beam-forming networks and apply matching 
to the decoupled ports of the antenna array. 

In this paper, we investigate optimal transceiver design 
for compact arrays from a communication theory perspec- 
tive. Specifically, we focus on the impact coupling has on 
the bandwidth of coupled circular arrays and derive optimal 
broadband matching networks. It will be shown that for 
broadband systems with uniform circular arrays, mutual cou- 
pling decomposes the coupled array with spectrally-identical 
spatial modes into spectrally non-identical eigen-modes - with 
differing bandwidths and resonant-frequencies. Similar results 
have appeared for specialized cases recently f T2l , lfT3l . These 
studies treat the problem in detail from a microwave theory 
standpoint, including many of the implementation aspects. Our 
findings and analysis help generalize how coupling impacts the 
RF bandwidth in the context of circular arrays. By combining 
concepts from antenna/microwave theory, Fano's broadband 
matching and Shannon's information theory, we present a 
unified communication-theoretic framework and evaluate the 
diversity limits of coupled broadband systems - utilizing 
orthogonal frequency division multiplexing (OFDM) - with 
varying antenna spacing. 

The organization of the paper is as follows. Sec. II presents 
an overview of the basic microwave theory tools necessary to 
analyze the problem of coupled MIMO systems. Sec. Ill illus- 
trates how coupling impacts the RF bandwidth of a compact 
circular array and introduces the concept of virtual antennas. It 
also discusses the applicability of Fano's broadband matching 
theory to characterize the optimal broadband matching net- 
work. Sec. IV supports these findings with numerical data. In 
Sec. V, we present a transceiver model for a broadband coupled 
array, and employ it in Sec. VI to develop a system model for 



0000-0000/00$00.00 © 2013 IEEE 



2 



JOURNAL ON SELECTED AREAS IN COMMUNICATIONS 



Diversity-OFDM in the presence of mutual coupling. Sec. VII 
presents results for the outage capacity as a function of antenna 
separation. We conclude by summarizing the main findings of 
the paper in Sec. VIII. 

II. Coupled Transceivers 

As the antennas in an array are brought closer, currents 
flowing in one element induce voltage across the other This is 
commonly referred to as mutual coupling. This phenomenon is 
usually modeled using the circuit (or network) representation 
of the antenna array. For an N element antenna array, the 
currents i flowing through the antenna ports, and voltages v 
induced across them can be modeled by TV x 1 vectors and 
N X N matrices as: 

V = Zai + Vo . 

Here, Za denotes the N x N antenna array impedance 
matrix and Vo the iV x 1 open-circuit voltage induced by 
the incident electro-magnetic field (EM). The diagonal entries 
of Za denote self-impedance and the non-diagonal ones 
mutual-impedance, such that in the absence of coupling, the 
antenna impedance matrix is diagonal. It has been shown 
that optimal performance can be achieved by transforming 
the coupled antenna array to an uncoupled one by use of 
a 2A^-port impedance-transforming network called matching 
network Zm, inserted between the array and the rest of 
the transceiver H, Q. The choice of the matching network 
depends on the antenna separation (essentially the extent of 
coupling) and is usually chosen as a fixed passive network. The 
optimal networks are generally non-diagonal in nature, while 
the more practical, but sub-optimal, are diagonal. However, 
these studies have essentially focused on narrowband models 
and in order to incorporate bandwidth considerations, the 
analysis must be extended to using matching networks that 
perform well over a given bandwidth. As will be shown 
later, the system analysis of broadband systems with mutual 
coupling is eased by the use of scattering-parameter or S- 
matrix representation, instead of the impedance matrix. 

To begin, we introduce the basic elements of broadband 
matching theory for a single antenna (or 2-port network) 
system. Consider a broadband 2-port network (shown in 
Fig.[T]i, with impedance matrix Z(s) terminated into reference 
impedances zi{s) and Z2{s) on either side. Here, s = a + juj 
is the Laplace variable and cj = 27r/ denotes frequency 
in radians/sec. For this simple (assuming reciprocal) 2-port 
network, the elements of S-matrix represent the reflection and 
transmission coefficients Fi, r2 and T, respectively, as shown 
in Fig. [T] 



S(.) = 



ri(s) 

T{s) 



T(s) 



It relates the voltages across and currents through the ports to 
the incident and reflected normalized wave vectors a and b, 
respectively, via lfT4l 
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Fig. I. Input and output wave vectors for a two-port network 



where. 



a(.) 



ai(s) 
02(3) 



b(.) 



bi{s) 
b2{s) 



For impedances = 22(5) = 1 il, the S-matrix is readily 

computed using the impedance matrix 



S(.) = (ZGs)+I)-i(Z(,s)-I) 



(2) 



where I denotes the identity matrix. For a lossless and recip- 
rocal network, S satisfies 

SO-c.)S^(jc.)-I, S(jc.) = S^O'c.) . 

In the next section, we discuss a very special but practical class 
of planar antenna arrays - uniform circular arrays (UCA). We 
shall use the circulant nature of a UCA to ease the analysis and 
gain rich insights into the impact of coupling on the bandwidth 
of compact arrays. 

III. Virtual Antennas 

We begin with characterization of the broadband antenna 
array S-matrix S^. In order to establish S^, it is convenient 
to look at the impedance matrix Za- For example, an antenna 
array with N — 2 identical elements placed a distance d apartQ 
has a symmetric impedance matrix of the form (cf. Fig. |2]i 

2l2(jw) Zii{juj) 

The symmetric nature of Za above enables us to express it in 
terms of its eigen-value decomposition (EVD) 

ZaUuj) = QAa{jlo)Q" , (3) 
where the set of unitary eigen-vectors is given by 



Q 



1 

71 



and the eigen-values by 

Kj'^) = {^ll(j^) + ^12(jw), Zii{juj) - Zi2{juj)} 



(4) 



h{s) = S(s)a(s) , 



(1) 



henceforth referred to as eigen-impedances. It is important 
to point out that the unitary transformation Q above is 
frequency-independent. This makes the analysis of the broad- 
band matching problem applicable to coupled arrays, and the 

'The antenna separation d is specitied in terms of Ac - wavelength 
con'esponding to the center frequency /c. 
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implementation of matching network considerably easier As 
will be discussed later, the matching network can be realized 
as a cascade of a frequency-independent decoupling network 
followed by a diagonal broadband matching network. 

A. Uniform Circular Arrays 

A careful observation of the structure of and the unitary 
transformation Q reveals that it is straightforward to extend the 
entire theory to uniform circular arrays. A uniform circular ar- 
ray has a circulant impedance matrix. For a circulant with 
N elements, the eigen-vectors for Z^(ja;) = QA^(ja;)Q^, 
are given by the columns of the unitary matrix 



1 a 



1 



N-l 



1 



a 

T,2(JV 



AN~1)(N-1) 



where a = e~'^'^^l^ . The eigen-values (A^) are given by the 
discrete Fourier transform (DFT) of the first row of T^a- 

For example, the spatial unitary transformation that decou- 
ples an 3 UCA, is given by 

1 1 

1 



1 



1 



1 

1 

"5 



V3 



1 , ,V3 
1 ,V3 



1 

and the eigen-impedances by 

Ai(jcj) = zx\{juS) + 2zi2(jw) 

hU^^) = A2(jw) . 
Similarly for A^ = 4 antennas, we have 

1 



(5a) 
(5b) 
(5c) 



Q 



1 1 

1 -J 

1 -1 

1 .] 



1 

-1 
1 

-1 



J 

-1 



and eigen-impedances as 

Ai(jw) = zii{juj) 
A2(jw) = Zii{juj) 
hU^^) = ziiijuj) 
A4(jw) = A2(jw) . 



2212 (jw) + Zl3{ju>) 
Zl3{j(^) , 

22:12 (jw) + Zi3{juj) 



(6a) 
(6b) 
(6c) 
(6d) 



The unitary transformation Q, is essentially an orthogonal 
beam-forming matrix implemented using RF networks called 
beam-formers [11 J. These beam-formers are capable of pro- 
ducing N spatially orthogonal beams, hence the operation 
represents a spatial DFT|3 In terms of the antenna radiation 
patterns, this operation can be thought of as decomposing 
the composite array pattern into overlapping, but mutually 
orthogonal patterns at the output of the beam-former. For 

"One of the well known and widely u.sed matrix for arrays with 2" number 
of antennas is the Butler matrix 1 15 |, representing spatial FFT operations. 
However, our case necessitates a more general decoupling matrix. 



^11 ^12 

^12 -^u 



^11 ~Zl2 



^12 



^11 ~ '12 



Fig. 2. Impedance matrix representation of a 2-element array 



an illustration of such pattern decomposition, see lfT2l . ifTSll . 
Similar analysis, however, can not be extended to uniform 
linear arrays for which the array impedance matrix is complex- 
symmetric 

Although the antenna array is an A^-port network, it can be 
appropriately extended to a 2A^-port network for mathematical 
convenience Q, such that 



S22a(jw) 
S21a(jw) 



S2lQ(jw) 
S22a(iw) 



where the original A^-port representation of the antenna array 
is represented by 8220 block, computed using 

S22a(.?^) = (ZaO'c^) + irHZAijio) - I) . 

The other blocks must be evaluated based on the lossless 
(S^S^ = I) and reciprocal (S^i = S^) properties of the 
antenna array. The symmetry of the system under considera- 
tion and that of the individual blocks constituting Sa, further 
allows us to write (using EVD) 



Q 

Q 







A21a(jw) 
A22a(jw) 











where 



A22a(j^) = (A^(j^)+I)-1(A^(JC.)-I) 



(7) 



B. Impedance Characterization 

Next, we characterize the individual entries of antenna 
impedance matrix Za over a broad range of frequencies. 
Fano's broadband matching theory requires an impedance to 
be analytic over the entire s-plane. In order words, the resistive 
part must be an even function of uj and the reactive part an 
odd function of lu. Since the eigen-impedances are a linear 
function of self and mutual-impedances, it suffices to consider 
an impedance model fit for the eigen-values of the antenna 
impedance matrix. 

Although the resistive part of the eigen-impedances is, in 
general, frequency-dependent, for reasons outlined next, we 
assume that it is fairly constant over the frequency range 
of interest (say, 10% relative bandwidth). First, the resistive 

^Although a unitary transformation for such matrices exists for each data- 
point over the frequency range, it is, in general, frequency-dependent. 
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Antennas 



Matching Network 



Load 




Fig. 3. Optimal matching implementation using decoupling networks 



part variation over frequency - as obtained from numerical 
EM code (NEC) simulations - will be shown to be much 
slower than the reactive part. Second, the system analysis 
is greatly simplified by such as assumption. The system 
transfer function, which is essentially governed by the antenna 
reflection coefficients (magnitude-squared), will be shown to 
closely approximate that of a series RLC modelQ Thus, eigen- 
impedances assume the form X{juj) — B. + jX{uj), where R 
and X represent the real and imaginary parts, respectively. 

The eigen-impedances can essentially be regarded as virtual 
antennas comprising an uncoupled antenna array. In order to 
analyze their frequency responses, we must first characterize 
the eigen-impedances. Consider an array of two ideal and 
identical A/2 resonant dipole antennas placed sufficiently far 
apart. Such a resonant antenna can be well modeled by a series 
RLC circuit within a 10% relative bandwidth EH, IHl. For 
such an array, is diagonal. A resonant antenna with a series 
RLC equivalent, can alternatively be expressed in terms of its 
resonant frequency luq and quality factor Q, where; 



1 



Wo 



, Q 



As the spacing between the antennas decreases, EM inter- 
actions start to alter the spectral responses of the individual 
antennas. Numerical simulations using NEC suggest that the 
eigen-impedances also conform to that of a resonant antenna. 
Hence, we model them using a series RLC equivalent, or 
resonance parameters {Q^loq) as 

A„(s) = Rn + LnS + — — , Rm in, C„ > , 
where L„ = RnQn/ojQn , Cn = l/QnRn(^On- The COITC' 

sponding reflection and transmission coefficients 

" r[{s) 



A„(ja;) = Rn 



1 + jQn 



Aiia(s) 
A21a(s) 



T',{s) 

Ti{s) 
ms) 



Without loss of generality, more rigorous impedance models can be 
employed to analyze such coupled MIMO systems using the approach outlined 
in this paper. 



normalized to i?„ fl impedances are given by 

1 + (s/a;o„)2 



1 + (,s/wo„)2 + 2(s/(3„wo„) 

-1 



(8) 



We refer to the frequency response of these virtual antennas, 
denoted by 

4/2 



as eigen-modes of the coupled array. 



(10) 



C. Broadband Matching Constraints 

The broadband matching network design for this uncoupled 
system is rather straight-forward. The proposed matching 
network implementation is illustrated in Fig. [3] The purpose of 
matching is essentially shaping these spectral responses such 
that the overall system has a frequency-flat response over the 
bandwidth of interest. Ideally, one would expect to choose a 
matching network that ensures 

|r„(/)|2 = i-|r„(/)|2, / e B , 

is unity. Here, Tn, r„ correspond to the transmission and 
reflection coefficient of the cascade of the antenna array and 
the matching network. However, Fano's broadband matching 
theory 11181 reveals that there exist gain-bandwidth trade- 
offs for physically realizable matching networks, built using 
lumped passive elements. It imposes certain integral bounds - 
determined by the source and load impedances connecting the 
matching network - on the matching efficiency of the network. 

For the series RLC model considered in our work, the set 
of broadband matching constraints are given by 

1 l^On \ ^ 



(a) / log 

IB 



|r„(/)p 



df 



(11a) 



(b) / /-'log 



|rn(/)|^ 



df 



Ev«. (lib) 



where Zr represent additional zeros in r„, that may sometimes 
be necessary to introduce in the right-half complex plane 
(Re(zr) > 0) in order to satisfy all of these constraints]! 

^ Zr must occur in conjugate pair if they are complex. The authors would 
like to point out a typo in 1191 . Eq. (14b); the correct bound is as stated above. 
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Normalized frequency (f/f ) 0.5 1 1 .5 

Normalized frequency (f/f ) 



Fig. 4. Impedance: N = 2, d = 0.25Ac 

Fig. 5. Antenna eigen-modes: N = 2, d = 0.25Ac 



Observe how the bound is inversely proportional to the antenna 
Q. Clearly, an antenna with a broader frequency response 
(lower Q) offers better gain-bandwidth trade-offs. 

IV. EiGEN-MoDE Illustrations 

Next we present numerical results for iV = 2, 3, 4 element 
array of dipole antennas with lengtlH 0.475Ac and radius 
for a variety of antenna spacings d. 

A. Impedance Parameters 

Fig. m shows the broadband eigen-impedances calculated 
using NEC for N — 2 and d — 0.25Ac. The resistive and 
reactive parts of the impedance are in general, non-linear 
functions of frequency. The figure also shows the imaginary 
part of eigen-impedances obtained by series RLC-fit which 
evidently, is in close agreement for bandwidths on the order 
of 10%. Although we consider a series RLC model for ease of 
analysis, the main results and ideas conveyed in this work can 
be applied to more rigorous models, as long as these models 
are analytic over the entire complex plane, in order to apply 
Fano's matching theory. 

Fig. |5] shows the corresponding eigen-modes found by 
curve-fitting the resonance parameters {Q,ujo) for the same 
setting. It clearly shows the impact coupling has on the 
bandwidth of the two virtual antennas. Mutual coupling is 
seen to decompose a two antenna coupled array with identical 
spatial modes into two spectrally non-identical eigen-modes - 
one broadband and the other narrowband. At much smaller 
spacings (e.g., d = 0.1 Ac), the narrower mode is essentially 
non-existent. The shrinking bandwidth of certain eigen-modes 
with increasing coupling is a manifestation of array super- 
directivity 1201 . Table H] summarizes this data for N = 2, 
d = 0.25Ac. 

Similar results have appeared in [12], where the common 
and difference modes (essentially, the eigen-modes) of a 2- 
antenna coupled array have been illustrated for d — 0.1 Ac, 

*The lengtli is chosen such that each antenna in isolation, has a relative 
resonance frequency of 1. 



along with the orthogonal radiation patterns and a simulation 
test-bed to realize the system. 1 13| further extends the findings 
and implementation results to a 4-antenna circular array. It also 
proposes equivalent circuit models for each of the four eigen- 
modes using ladder LC networks, as opposed to the series 
RLC fit assumed in this paper. 

B. Usable Bandwidths 

Microwave/RF bandwidths are usually parameterized by 
voltage standing-wave ratio (VSWR): a metric based on 
matching efficiency, indicative of the range of voltage fluctu- 
ations in the standing wave formed due to reflections arising 
from an impedance-mismatch ||2TI . 



VSWR 



1 + \n^) 

i-irH 



The higher this ratio, the larger the mismatch and smaller the 
bandwidth. Note that < |r(a;)| < 1, implies VSWR > 1. 
A convenient measure of the RF bandwidth is defined as the 



TABLE I 

ElGEN-IMPEDANCE PARAMETERS: : 



2,d = 0.25Ac 



Parameter 


Value 


Antenna type 


Dipole 


Antenna length 


0.475Ac 


Antenna radius 


O.OOlAc 


Quality factor (Qi) 


3.75 


Resonant frequency (/oi) 


1.0425/c 


Ri 


118.76 n 


Li = QiRi/ujai 


67.99//c H 


Ci = l/QiRiujm 


342.78//c fjF 


Quality factor {Q2) 


16 


Resonant frequency (/02) 


0.9675/c 




28.31 n 


L2 = Q2R2/^02 


74.53//c H 


C2 = 1/Q2-R2l^02 


1.4//c mF 
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frequency range such that |I2T1, US), 1 < VSWR < 2, i.e., 

|r(w)| < 1/3 =^ 1 - |r(w)p > 8/9 . 

Fig. |6] shows the VSWR and usable bandwidth for eigen- 
impedances calculated using NEC with termination Ri and 
i?2, respectively. 

C. Discussion 

These results shed substantial light on the impact coupling 
has in governing the RF bandwidth of compact arrays. It 
offers intuitively-coherent interpretations - as the spacing goes 
to zero, in limit, the narrowband mode must vanish, leaving 
behind a system that behaves like a single antenna system. 
The other key observation - that the two modes are centered 
at different resonant frequencies - explains that a multiport 
matching network designed with narrowband assumptions is 
sub-optimal and that physically realizable optimal networks 
ought to be optimized over the bandwidth. It also highlights the 
importance of antenna design that can exploit this phenomenon 
by trimming the antenna lengths and centering these modes so 
as to maximize the diversity gains. 

Fig. |7]-[9] illustrate eigen-mode behaviors for different uni- 
form circular array sizes N = 2,3,4, and antenna separa- 
tions As expected, the impact of coupling at smaller spacings 
becomes profound as the array size grows. An analogy for 
the asymmetric behavior of eigen-modes can be drawn from 
quarter-wave transformers, a well known matching technique 
in microwave literature [22]. Each virtual antenna can be 
thought of as a transformer with a different length and 
impedance, thereby matching the other virtual antenna at a 
different frequency with a different bandwidth efficiency. 

V. Diversity-OFDM for Compact Arrays 

Having analyzed how coupling impacts the RF bandwidth 
in the context of circular arrays, next, we evaluate the capacity 
of broadband diversity systems for variable antenna spacing. 
We choose OFDM as the broadband transmission scheme for 
our system. But first, we introduce the basic elements of 
multiport networks useful in modeling the coupled channel 
vector, followed by a system model. 

A. Broadband Coupled Array Model 

To that end, consider the receive diversity system with 
N antennas, as shown in Fig. [TO] in its S-matrix network 
representation. It shows the cascade of two 2Af-port networks 
- Na, representing a coupled lossless and reciprocal antenna 
array , and Nm representing the lossless and reciprocal 
matching network - terminated into a bank of uncoupled 
load impedances z^. The load here is indicative of low 
noise amplifiers and other downstream components of an RF 
chain, primarily, mixers and A/D converters. The EM field 
incident on the receive antenna array induces an open-circuit 
voltage across the antenna terminals, which acts as the source 
excitation, represented by the input wave vector ai toward the 
left of the antenna array. 

'in the context of a UCA, d is defined as tlie separation between adjacent 
antennas. 




1 1.05 1.1 1.15 

Normalized frequency (f/f ) 



Fig. 6. VSWR of a 2-eiement array for d = 0.25Ac 



The 2N X 2N S-matrices for the antenna array and matching 
network (in x block-matrix format) normalized with 
respect to 1 reference impedances are given by 



Slla Si2o 
S21a S22a 



Sllr 
S21r 



Sl2m 
S22m 



where we have omitted the frequency-dependence by sup- 
pressing (s) for aesthetic reasons. We shall henceforth assume 
it implied, unless stated otherwise. The cascaded 2A^-port 
network has an S-matrix 



Sllc 

S21C 



Sl2c 
S22C 



where, 

Sllc 
Sl2c 
S2IC 
S22C 



represents the cascading operation and 

Slla + Si2a(I — SiimS22a) ^SiimS21a 
Si2a(I — SiimS22a) "'"Si2m , 
S21m(I ~ S22aSllm) ^§210, 



(12a) 
(12b) 
(12c) 



= S 



22m 



S21rri(I — S22aSllm.) S22aSl2m 

.(12d) 



The inward and outward traveling wave vectors are related by 
([Til, except that, the input and output wave-vectors at the left 
are now vectors, denoted by ai and bi, respectively. Those 
on the right side are denoted by a2 and b2, such that 



bi 
b2 



Sllc Si2c 
S2IC S22C 



ai 
a2 



(13) 



It can be easily shown that optimal matching combined with 
the decoupling network (spatial-DFT) can be realized by 



Q 
I 



Aiim 
^21m 



Ai2m 
^22m 







" 







I 



such that from (fT2] i. the overall cascaded network can be 
rewritten in block-matrix format as 

Q ] r r T ] r 

Q T r 
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(a) d = O.lXc (b) d = 0.5Xc 



Fig. 7. Antenna eigen-modes: N = 2 





Fig. 9. Antenna eigen-modes: N = 4 
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Source/EM Excitation 



Cascaded Network 



Lx)ad 




Antenna 
Array (A',, ) 



'llii '-'12a 



'2la "220 



Matching 
Network (A',,,) 



S,2 



S2.,„ S 



12m 
22 ni 



b2(^) 



Fig. 10. S-matrix representation of RF front-end 



where, 

r = Alia 



Ai2a(I — Aii„iA22a) "'"AiimA21o 



From a broadband matching perspective, it suffices to analyze 
the constraints on the eigen-values of Suc given above. This is 
once again equivalent to saying that the problem of matching 
the coupled identical antenna aiTay Za has been replaced by 
that of a uncoupled non-identical virtual antenna array Aa- 
The entries of 

' ri(s) 

T^is) 

represent the reflection coefficients at the output of the match- 
ing network. 



r(.) = 



B. Signal Model 

We begin by presenting the signal model of a traditional 
Diversity-OFDM system with 1 transmit and N receive an- 
tennas, which we shall later extend to incorporate coupling. It 
is assumed that the separation between the antennas is such 
that coupling between them is negligible. 

It is well known that the use of orthogonal sub-caniers 
in OFDM with a cyclic prefix converts a frequency-selective 
MIMO channel into a set of paraUel frequency-flat MIMO 
channels ll23l . A Diversity-OFDM system with K sub-carriers 
(spanning bandwidth B) modulated by symbols Sk where k 
represents the fc-th sub-carrier, is well modeled by 



Tk = hfcSfc + nfc , fc = 1, . 



(14) 



Here, is the x 1 received vector symbol on the fc-th 
sub-carrier and is the A x 1 additive white Gaussian noise 
(AWGN) vector at the receiver with zero mean and covarianc^ 
R„^ = E[nfcnf ], denoted by nfe - C7V(0, R„^ ). The channel 
vector for the fc-th sub-carrier is given by 

hi[k] " 



fc = 1, 



.,K. 



^E[-] represents the expectation operator. 



The transmit and receive spatial fading-correlation is modeled 
using the Kronecker model |24| such that the Z-th tap time- 
domain channel vector (obtained via inverse-Fourier trans- 
form) can be expressed as 



(15) 



where h^,; represents the 1 x A white channel vector having 
i.i.d. coniplex Gaussian entries with zero rnean and unit 
variance h^; ~ CA'(0, 1), and Rh = (l/A)E[h(hP] denotes 
the receiver correlation. 

Owing to the orthogonal decomposition of the frequency 
selective channel, the cumulative Diversity-OFDM system can 
be represented in matrix notation by 

r = 'Hs + n, (16) 

where, is a KN x K block diagonal matrix given by 

hi 



n 



The KN X 1 received vector r. A" x 1 transmit vector s and 
KN X 1 noise vector n are given by 





ri 




Sl 




ni 


r — 




, s = 




, n = 





The Shannon capacity for such a system (in nats/s/Hz) is given 

by m 



C 



^ umx {logdct + n-^nits-H"^ } (17) 



where, / is an NK x NK identity matrix, and Ti-s and Tl-r. 
are the transmit-signal and noise covariances: 



(18) 



There have been numerous extensions to the above model in 
order to account for some of the channel non-idealities, such 
as fading correlation at the transmitter and/or the receiver, 
either due to smaller antenna separation or non-richness of 
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the multipath fading environment. The broadband antenna 
behavior, however, has mostly been assumed ideal. In the 
next section, we introduce a channel model that incorporates 
antenna coupling into the signal-model outlined in ( fT4b . 

C. Coupled Signal Model 

To model the impact of coupling, hk must be modified 
to include the transmission matrix or transmissivity of the 
cascaded network (of antenna array and matching network) 
modeled by S2ic(/fc). We skip the tedious network analysis 
and present the effective channel vector at the fc-th sub-carrier 
directl)0: 

= S2ic.fchfe . (19) 

In the presence of mutual coupling, the signal model can thus 
be expressed for the fc-th sub-carrier as 



rfc = SfchfeSfe + rifc , 
where Sfc = S2ic,fc = S2ic(/fc)- 



(20) 



where, ks is Boltzmann constant, B is the system bandwidth, 
= Re(Zyi) is frequency independent (by assumption), 
and Tc represents the correlation between forward and reverse 
traveling noise. Due to the limited scope of this paper, we 
assume Tc is negligible compared to Tf and T^, and restrict 
the impact of various noise sources by varying Ta, Tf and Tr 
relative to the standard temperature Tq. 

Observe that Sfc and Ra admit the familiar EVD 

Sfc = QTfcQ^ , Ra = QRc(A^)Q^ . (22) 

This allows us to diagonaUze the noise covariance by Q, i.e., 

Rrifc = QSn^Q^ , 

such that, 

S„, = AkBB[{TA~Tr)'R.e{AA){l-TkT^) + ... 



(23) 



where we have used the lossless property of the network, i.e., 
TkT^k '^k^k ■ normalize the noise covariance such 



that for i.i.d. case, 



iVo = 4fcBS [Ta Ke{zA){l ~ \Tud?) + Tf + T^lr^.^p) 



D. Receiver Noise Model 

The additive noise nfc is usually modeled as a combination 
of noise from various sources in the RF chain Q- In general, 
the noise sources can be categorized into three types: (a) sky 
noise or antenna noise, consisting of thermal radiation, cosmic 
background, and interference from other devices, (b) amplifier 
noise, and (c) downstream noise, consisting of noise from the 
rest of the RF chain components. 

Alternatively, we classify the noise as antenna noise, and 
load noise - a combination of amplifier and downstream noise 
IIZTI . Furthermore, the load noise in general can be considered 
to be a combination of forward traveling noise n/, and reverse 
traveling noise n^. Thus, the total noise at fc-th sub-carrier 
referenced to the load, is given by 



nfc = S2ic,fcns,fc + n/^fc + S22c,fcnr,fc . 



(21) 



The sky noise and load noise can be well modeled as statisti- 
cally independent, zero-mean, circularly symmetric, complex 
Gaussian (ZMCSCG) and spectrally white: 

Us ^ CU{0,'ikBTAB'R.A) , n/ - CA/'(0,4fcBT/BI) , 
^CJ^{0,4:kBTrBI) , 

a reasonable assumption for bandwidths less than 10%. Here, 
Ta denotes the antenna temperature in Kelvin, while Tf and 
Tr are the effective noise temperatures which can be computed 
from the amplifier noise parameters (cf. |28, Chap. 1]). 

The reverse and forward traveling noise waves are in 
general, correlated to an extent determined by exact amplifier 
models (cf. Il28l Chap. 1]). and the noise covariance by 



R„. = 4kBB[TASkRASj^ +TfI+... 

... + r,(I-SfeSf)-2Re(T;Sfc)] 

'interested reader is referred to 1261 . 



E. Capacity 

The cumulative Diversity-OFDM system in the presence of 
mutual coupling can be written in matrix format (similar to 
(HSll), as 



r = STis + n 



(24) 



where S is the KN x KN block-diagonal matrix given by 

" Si 



The Shannon capacity that incorporates the impact of mutual 
coupling can thus be represented by 

C = 4 max {log dct ( / + ^U:^^ SHUsU" S" 
K -R-^,s [ ^ V ^ 

where the optimization space now also includes the matching 
network design. A recent study has analyzed information- 
theoretic limits of such a system in the presence of channel 
state information (CSl) that jointly optimizes transmit power 
allocation and receiver broadband matching |29 | - such that 
the optimal solution follows a mutual space-frequency water- 
pouring characteristic. 

In this work, we address a receiver design that operates 
fairly well over the entire signal bandwidth and is indepen- 
dent of the power allocation and channel fading conditions. 
By employing appropriate decoupling networks, the capacity 
subject to uniform power allocation across K sub-carriers 
(TZs = Sgl), can be expressed as 

1 / £ 

C = — maxlogdet / + -f n;^^S'H'H"s" 

K s \ iVo 

= ^ ui&yi log Aei ( I + ^ n" QY.-^ TT" Q^U 
K s \ Nq , 
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where we have expressed S and TZn as: 

5 = QTQ" , TZr^ = QS„S^ , 
using KN X KN block-diagonal matrices 



Q = 



Q 



Q 



Ti 



Tk 



The capacity can thus be simplified to 



N, 







r;crf)h. 



where, = Q^h^ represents the effective fading path-gains. 

A simple and practical solution to the problem at hand is a 
box-car matching characteristic at the receiver, defined as; 



r(/) = 



To, 
1, 



elsewhere 



such that 

C = 



1 ^ 

k=l 



log i + -f hfs„i(i-rfcrf)h. 



VI. Results 

Monte-Carlo simulations are carried out for 100, 000 chan- 
nel realizations. We assume a quasi-static (or block) fad- 
ing channel, i.e., the channel remains constant during each 
OFDM symbol. The fading paths gains are modeled as i.i.d. 
complex Gaussian entries with zero mean and unit variance, 
hwi ^ CJ\f{0, 1). For an iV-antenna uniform circular array, the 
incident electric field is modeled in NEC as a superposition 
of K' = 32 vertically polarized plane waves with i.i.d. phases 
uniformly distributed on [0, 27r). The angles-of-arrival (AOA) 
of the plane waves, (f>Q, . . . , (j)K'^i, are uniformly spaced in 
azimuth from to 27r. Under these conditions, the open-circuit 
fading path gains for m-th and n-th antenna separated by d„in 
are approximately Gaussian with correlation matrix (cf. (fTSb ) 



[Rh] 



K'~l 

E 

k=0 



where gn{4') is the open-circuit voltage induced in the nth 
antenna by a zero-phase plane wave with AOA normalized 
so that J2k \9n{4'k)\'^ — 1 for an isolated dipole. 

For each virtual antenna Tk, we consider a box-car match- 
ing profile over a relative bandwidth W — B/fc, such that 
the matching constraints ( fTTT i manifest themselves as 

(a) W^log ^ ^ ^ _ 1 z„ , (25a) 



ib) 



iroP 
w 



log 



(1-1^2/4) "irol 



Q fo^ 

1 27r 



1 

To 



(25b) 



-B — With Matching 
-e — Without Matching 




0.1 



0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 
Antenna spacing (d/A. ) 



Fig. 11. Diversity-OFDM: Outage capacity vs. spacing, N = 2 



It suffices to consider a pair of complex zeros in the right-half 
complex plane z^i = z*2 iflTjl (see Appendix), such tha0 



Wlog 



ro 



< (i - W^/A) 



27r 



Q 



and 



-27r{l-W^ /4)/QW 



The usable eigen-modes are determined based on the VSWR 
< 2 criteria. By adjusting the antenna lengths, the two eigen- 
mode resonant frequencies can be altered to maximize the 
overlap around /c, so as to increase the diversity order of the 
system. The OFDM parameters for simulation are inspired by 
IEEE 802.11a standard, with K = 6A sub-carriers spanning a 
bandwidth of i? = 20 MHz. The results are presented for a 
load-noise dominant scenario: {Ta ■ Tf : Tr) = (1:2:0), for 
W = 2% and Ss/Nq = 10 dB SNR. We measure the system 
performance in terms of the outage capacity - probability that 
the capacity falls below a certain threshold Cq: 

Cout = Pr(C < Co) . 

Fig. [TT| shows 1% outage capacity for different antenna 
spacings for two element arrays. The performance predicted by 
these simulations shows that with increasing antenna spacing, 
the capacity increases for small spacings and starts to saturate 
around a quarter of a wave-length spacing. This behavior 
is in stark contrast to some of the studies which (assuming 
narrowband models) predict that, in theory, it is possible to 
achieve IID performance using optimal multiport matching 
networks even at very small spacings |7|. The answer clearly 
lies in bringing the bandwidth variable into the equation 
and re-visiting the communication-theoretic formulation for 
MIMO systems with mutual coupling. The results also suggest 
that in order to exploit diversity, practical arrays can be built 
with a separation of less than half-a-wavelength spacing - an 
oft-cited assumption in MIMO literature. 

'"Strictly speaking, considering only the first constraint yields an upper 
bound on capacity via IFoP = e"^'^/'^^, while only the second constraint 
yields a lower bound via IToP = e-2'^(l-'^^/4)/'3W'. For small W, the 
upper and lower bounds are quite tight. 
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Fig. [12] shows 1% outage capacity vs. antenna spacing for 
a four element UCA at Ss/Nq = 10 dB SNR, for the same 
settings. The impact of relative noise level from various noise 
sources, can be studied in a similar manner by varying Ta, Tj, 
and Tr relative to To; the capacity still grows with increasing 
antenna separation and eventually saturates to i.i.d. case. 



VII. Conclusion 

In this paper, we primarily investigated the relationship 
between coupling and the bandwidth of compact multi-antenna 
systems, and derived the capacity limits of a diversity system 
with box-car matching. We analyzed the spatial modes of a 
compact antenna array and demonstrated that they exhibit dif- 
ferent bandwidths and resonant frequencies; a phenomenon not 
observed in arrays with large antenna separation. Furthermore, 
we showed that traditional broadband matching theory can be 
applied to these eigen-modes in a straight-forward manner for 
uniform circular arrays. Although we considered equal band- 
widths for all of the eigen-modes, a differentiated approach 
with respect to the mode bandwidths can be applied to further 
optimize the performance of Diversity-OFDM systems. 

We discussed, how, in limit as the spacing between the 
antenna elements goes to zero, the most narrowband mode 
vanishes under strong coupling, leaving behind a system with 
a lower diversity order The other key observation that the 
these modes are centered at different resonant frequencies 
- explains that a multiport matching network designed with 
narrowband assumptions is sub-optimal and that physically 
realizable optimal networks ought to be optimized over the 
bandwidth. 

We also presented a communication-theoretic framework for 
Diversity-OFDM systems with mutual coupling and broadband 
matching. The results show that capacity increases with an- 
tenna spacing and that a quarter of a wavelength separation 
might suffice for most practical applications. An information- 
theoretic approach unifying the transceiver design with the 
overall system design, including signal processing aspects can 
lead to a new theory of compact MIMO communications. 



Appendix 
Matching Constraints 

Let us consider a pair of complex zeros in the right-half 
complex plane z^i — z*2 = a + jp. Substituting normalized 
frequency /„ — f / Jq in (fTTI) . the matching constraints are 
given by 



1 



TOP 
1 



df 



^ |r(/)p 



df 



Q 

2nfo 



{Zrl + Zr2) 



Q 

27r 2a 



2a 




0.4 0.5 0.6 0.7 
Antenna spacing (d/A.^) 

Fig. 12. Diversity-OFDM: Outage capacity vs. spacing, N = A 



and 



log TFTKH df = - {Zrl + Zr2) 



T2 log 



\nf)? 

1 



|r(/)P '^^ foQ \Zrl\' 



loqQ 



2a 



2tt 2a/o 
dfn = -pr 



/n^^°^|r(/„)P Q \Zrir 

For a box-car reflection coefficient, r(/„) — Tq over the 
frequency-range 

l-W/2< fn<l + W/2 , 
the two constraints yield: 



WGq = 



27r 2a 



W \ ^ 27T 2afo 

where Go = — log jFop. By choosing jz^ip = /q, we get the 
least upper bound 

W 



1 - J 



Q fo 



Choosing a >• /3 (such that a ~ vTo and a/fo<S^ 1), we 
have 



Go < (1 - W^/i) 



2tt 
QW 
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